LIFSHITZ TAILS FOR MATRIX- VALUED ANDERSON MODELS 



H. BOUMAZA AND H. NAJAR 

Abstract. This paper is devoted to the study of Lifshitz tails for a continuous matrix- valued Anderson- type 
model i/u acting on L^(K'') ® C^, for arbitrary d > 1 and D > 1. We prove that the integrated density of 
states of Hu, has a Lifshitz behavior at the bottom of the spectrum. We obtain a Lifshitz exponent equal to 
— d/2 and this exponent is independent of D. It shows that the behaviour of the integrated density of states 
^ ' at the bottom of the spectrum of a quasi-d-dimensional Anderson model is the same as its behaviour for a 

^ ^ , d-dimensional Anderson model. 

> 

o 

1. Introduction 

1.1. A general model. We study the Lifshitz tails behaviour of the integrated density of states (IDS for 

(— I . short) of random Schrodinger operators of the form: 
» 



> 



(1.1) Ho{co) = -Ad ® Id + J] yj")(x - n), 

acting on the Hilbert space L2(M'^) C°, where d > 1 and D > 1 are integers, A^^ is the d-dimensional 
Laplacian and Id is the identity matrix of order D. 

Let (O, a, P) be a complete probability space and let oo ^ Q. We assume that, for every n G Z"^, the functions 
X I— >■ vS'^\x) take values in the space Sd(IK) of real symmetric matrices of order D and that these functions 



■ are supported on [— i, and bounded uniformly on x, n and uj. The sequence (Kj" )„g2'' is a sequence of 
^ " independent and identically distributed {i.i.d. for short) random variables on Vl. We finally assume that the 
sequence (Kj"^)„gzd is such that the family of random operators {iIo('^)}wef^ is Z^-ergodic. An operator 
like (jl.ip is also called a quasi-d-dimensional Anderson model. 
The vector space L^(R°') (8) is endowed with the usual scalar product: 

< /,5 >L2(Rd)$5CD= / < f{x),g{x) >cD dx = V / fi{x)gi{x)dx, 

where / = (/i, . . . , /d), g = {gi,. . . , gn) and fi G L^(R'^) 'S>C, gi e L'^{W^) (g) C are the i-th components of / 
and g. 
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As the functions x vS'^\x), for every u ^ ft and every n € Z*^, take values in Sd(K) , the operator Hq{ui) 
is self-adjoint on the Sobolev space i?^(R°') ^C^, for every cj G J7. Thus, its spectrum a{Ho{uj)) is included 
in M. Moreover, due to the hypothesis of Z'^-ergodicity of the family {ifo(w)}tjGn) there exists a set Sq C M 
with the property: for P-almost-every u G fl, a{HQ{uj)) = Sq. 

1.2. Existence of the IDS. We want to study the asymptotic behaviour of the IDS associated to Hq{ljj) 
near the bottom of the almost-sure spectrum of Hq(uj). The IDS of Hq{uj) is the repartition function of 
energy levels, per unit volume, of Hq{uj). To define it properly, we first need to restrict the operator Hq{oj) 
to boxes of finite volume. We set, for L > 1 an integer, 

d 



(1.2) Cl 



2L + 1 2L + 1 



Then, we consider Hq^Cl{'-^) the restriction of Hq(uj) to the Hilbert space L^{Cl) O C° with Dirichlet 
boundary conditions on the border dCi of Cl- To define the IDS, we now consider, for every € M, the 
following thermodynamical limit: 

(1-3) No{E) = ^hm^ J^lTW* {X < E \ X e a{H^,c, (w))} • 

We have already proved in [Tj that, for the general model Hq{uj), for every G M, the limit ()1.3p exists 
and is P-almost-surely independent of a; G (see [H Corollary 1]). The question of the existence of (|1.3p 
involves two problems to solve. First we had to prove that, for every G M and every cj G ^2, the cardinal 
#{A < £^1 A G a{H()^CLi^))} is finite. Then we had to prove the existence of the limit when L tends to 
infinity. Both solutions to these two problems rely strongly on the fact that the semigroup (e~*^''''-^i'^'^'')f>o 
has an L'^-kernel, which is given through a matrix- valued Feynman-Kac fomula (see [1, Proposition 1]). Once 
we obtain that the cardinal #{A < -E| A G (t(//o,Cl ('^))} is finite, we prove the convergence, as L tends to 
infinity, of the sequence of Laplace transforms of the counting measures of the spectral values of Hq^Ql i^) 
smaller than E. We prove the convergence of this sequence by using Birkhoff 's ergodic theorem which leads 
to the existence of a Borel measure no on M, independent of uj, which is the desired limit. We finally set 

(1.4) yE G M, No{E) = no((-oo, E]), 

the distribution function of riQ. The measure no is called the density of states of Hq{uj). 

1.3. A particular model. After this review of existence result of the IDS for the general model H(){uj), 
we may consider a particular example of such model for which we will be able to prove precise results on 
Lifshitz tails of the IDS at the bottom of the spectrum. 



We consider 
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(1.5) = -Ad + W{x) + J2 

acting on L'^{W^) (g) C^, where d > 1 and D > 1 are integers, A^; and Id are as in We set: 

(1.6) H = -Ad0lj) + W and V^x) = 

For the model (jl.Sp . we make the assumptions: 
(HI) W : ^ Sd(M) is Z'^-periodic, measurable and bounded. 

(H2) Vi, . . . ,V£, are nonnegative, bounded, measurable, real- valued functions supported on [— i, ^]'^. More- 
over, we assume that, for every i G {1, . . . , D}, there exists a non-empty cube Cj C [— ^, ^]'^ which is 
not reduced to a single point such that > 1^. , where 1^. is the characteristic function of Cj. 

(H3) For every i € {1, . . . ,D}, (<^j-"^)„g2d is a family of i.i.d. random variables on a complete probability 
space {Qi,Ai, Pi), which are bounded, and whose support of their common law contains zero and 
is not reduced to this single point. Moreover, we assume that, 

n '7^ w^n r>i r log | log Pi(tjf ^ < e)l 

(1.7) G {1, . . . , D|, limsup = 0. 

e-^o+ loge 

(n) 

In particular, we can take Bernoulli random variables for the cj) s. By adding a suitable constant diagonal 

(n) 

matrix to the periodic background W, we may always assume that the uj^ s are nonnegatively- valued 
(because of their boundedness) . If we set 

{n,A,P) = (g)(J7i®...0j^D), (g)(Il®---®lD), (g)(Pl®---^PD) , 

then {^},A, P) is a complete probability space and {-ffajjojeo is Z'^-ergodic because of the non-overlapping of 

(n) 

the random variables . We denote by S the almost-sure spectrum of Hi_j. By adding a suitable scalar 
matrix A/d to the periodic potential W, we may always assume that inf S = 0, where S is the almost sure 
spectrum of the Z'^-ergodic family {H^^y^^^Q. 

The model (11. 5p is a particular case of (II. ip for which the potential split into a deterministic periodic part 
W and a random part V^j which appears as a diagonal matrix. We will denote by : E ^ -^(-^) the IDS 

Remark 1.1. If we assume that, at least for one x G W^, W{x) is not a diagonal matrix, then we cannot 
write as a direct sum ^iLi-f^^.i of scalar-valued operators H^^i acting on L'^iW^) (S> C and for which all 
the results we will present here are already known. 
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Remark 1.2. The hypothesis {H2) on the boundedness of the functions Vi and the houndedness of their 
support implies in particular that each Vi is in L^(M'*) (8) with p = 2 if d < 3, p > 2 if d = 4 and p > j if 
d > 5. These are the assumptions made in |10] . 

For d = 1, matrix- valued operators as (jl.Sp are also called quasi-one-dimensional Anderson models. Local- 
ization results in both dynamical and spectral senses for such models, for particular simple choices of W and 
Vi, . . . ,V£), are obtained in [2] and [3j. These quasi-one-dimensional models are of physical interest as they 
can be considered as partially discrete approximations of Anderson models on a two-dimensional continuous 
strip. Such a two-dimensional Anderson model on a continuous strip can, by example, modelize electronic 
transport in nanotubes. Indeed, a two dimensional continuous Anderson model is defined by: 

(1.8) i/es(c^) = -A2 + J^c^("V(x-n,y), 

acting on L^(M x [0, 1]) ^ C with Dirichlet boundary conditions on M x {0} and M x {1}. The a;^"^'s are 
i.i.d. random variables and V is supported in [0, 1]^. As the continuous strip M x [0, 1] has one finite length 
dimension ([0, 1]) and one infinite length dimension (R), we can physically consider this strip as a quasi-one 
dimensional nanotube. The spectral properties of Hcsiuj) describe properties of the electronic transport in 
the nanotube M x [0, 1]. 

1.4. The behaviour of the IDS. The main result of this paper is about Lifshitz tails for the IDS N{E) of 
Hi_j at the bottom of the spectrum. In 1963, Lifshitz (see ^12]) had conjecture that, for a continuous random 
Schrodinger operator of IDS N{E), there exist ci,C2 > such that N{E) satisfies the asymptotic: 

(1.9) N{E) ~ ci exp(-C2(^ - ^o)"^), 

as E tends to Eq, where Eq is the bottom of the spectrum of the considered Schrodinger operator. The 
behaviour (11. 9p is known as Lifshitz tails (for more details, see part IV. 9. A of [18j) and the exponent —d/2 is 
called the Lifshitz exponent of the operator. The principal results known on Lifshitz tails are mainly shown 
for Schrodinger operators, in both continuous and discrete cases (see [Bl El [lOl [HI [171 [H] and others) and 
for Schrodinger operators with magnetic fields (see [HI [Hj ). Up to our knowledge, all studied examples 
of Schrodinger operators are for scalar-valued operators, and no adaptation of the known results to matrix- 
valued operators like (II. 5j) has been done yet. 

In a previous article of one of the author (see [I]), we had already obtain a result of Holder continuity of 
the IDS for a particular example of model H:_j, in dimension d = 1. For d = 1, we can use the formalism of 
transfer matrices and define Lyapunov exponents for Hi^j and, in this case, the sum of the positive Lyapunov 
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exponents is harmonically conjugated to the IDS of H^, through a so-called Thouless formula (see [H 
Theorem 3]). It allows us to prove that under some assumptions on the group generated by the transfer 
matrices, the Fiirstenberg group H^, we have positivity of the Lyapunov exponents, their Holder continuity 
with respect to the energy parameter and thus, by the Thouless formula, the same Holder regularity for 
the IDS. This assumptions are hard to verify for a general D (where D is the size of the matrix- valued 
potential), but we were able to verify them for a very particular example of H^^ in dimension d = 1, were 
the periodic and random potentials W and acts like constant functions. In [3j, we studied the following 
Anderson operator: 



(1.10) H,{u) = -^^I^ + V + Y^ 



acting on L2(M) (g, C°, where D > 1 is an integer. Id is the identity matrix of order D and ^ > is a 
real number. The matrix F is a real D x D symmetric matrix. The constants ci, . . . ,cd are non-zero real 
numbers. For / C M, 1/ is the characteristic function of /. The random variables w-""* are like in model (jl.Sp . 



As we can see, Hi{uj) is a particular example of H^^ with W constant and Vi = l[o,^]) for every «€{!,..., D}. 
For this operator we had obtain the following regularity result: 

Proposition 1.3. [3, Proposition 6.2] For Lehesgue- almost every V € Sd(IR) , there exist a finite set Sy C M 
and a real number £c '■= £c(D,V) > such that, for every i € {0,ic), there exists a compact interval 
I{T>,V,i) C M such that, if I C I{D,V,i) \Sv is an open interval, then the integrated density of states of 
Hi{uj), E I— 7- Ni[E), is Holder continuous on I. 

Remark 1.4. We actually proved even more : in such an open interval I with Sn / 7^ 0, we have Anderson 
localization in both spectral and dynamical senses. 

The Proposition 11.31 is interesting in itself but doesn't give any information about the behaviour of the IDS 
at the bottom of spectrum and, until now, it was not clearly stated that it has a Lifshitz behaviour. One of 
the motivation of the present article is to fill this lack of information of the IDS for quasi-one-dimensional 
operators and in particular those like Hi{ijj) we have studied before from the localization point of view. 

1.5. The result. We can now state the main result of the present article. 

Theorem 1.5. Let be the operator defined by 1 j.5|) and let N be its IDS. We assume hypothesis {HI), 
{H2) and {H3) and we also assume that inf S = 0. Then, 

In particular, this limit does not depend on D. 
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Remark 1.6. (1) Under some assumption on the behavior of the integrated density of states of the hack- 
ground operator H , it might he possihle to ohtain a result for internal hands. 

(2) Theorem \1.5[ could he used to give a different proof of localization than the one provided in ^ (see 

[a da [Ml [22];. 

It is important to insist on the fact that the Lifshitz exponent —d/2 obtained here does not depend on 
the integer D > \. It means that, looking only at the Lifshitz behaviour of the IDS at the bottom of the 
spectrum of the considered operator, we cannot distinguished a quasi-d-dimensional Anderson model like 
(jl.Sp from a d-dimensional Anderson model (for D = 1). 

One of the motivations in considering matrix-valued Anderson models is that we could expect, as D tends 
to infinity, that we could obtain informations about a (d + l)-dimensional Anderson model from a quasi- 
d-dimensional Anderson model. In particular, by obtaining a localization result for (|1.10p for an arbitrary 
Z) > 1, we could have expect to obtain a similar localization result for the continuous strip (jl.Sp . The 
presence of the Lifshitz tails behaviour of the IDS is usually a strong sign of the presence of localization at 
the bottom of the spectrum. If we wanted to use the Lifshitz tails behaviour of the IDS to prove localization 
(like in ^) and at the same time following the idea of approaching a, [d + l)-dimensional Anderson model 
by a quasi-d-dimensional Anderson model, we would have expected a Lifshitz exponent depending on D 
in a way such that this exponent would tend to —{d + l)/2 as D tends to infinity. But, Theorem 11.51 
contradicts this. So, we actually obtained an argument in favor of the idea that we cannot really get a 
proof of localization in dimension 2 (or more generally in dimension d+\) by an approximation procedure 
using quasi-one-dimensional Anderson models. This, at least if we follow a localization proof based upon 
the Lifshitz tails behaviour of the IDS. 



2. Matrix- VALUED Floquet theory 

2.1. Matrix-valued Floquet decomposition. In this section, we review the main results of the Floquet 
theory for the deterministic operator 

(2.12) H = -Ad ^Id + W 

acting on L'^{W^) (8i C^, with periodic potential W, and we adapt them to the matrix- valued setting. More 
precisely, we assume here that is a Z^'-periodic function in LP(M^) C° with p = 2 if d < 3, p > 2 if 
d = 4 and p > | if d > 5. If is Z'^-periodic, measurable and bounded as in model (jl.Sp . it is in such an 
L^(M'^) (8) space. Then, H is essentially selfadjoint on C^(R'^) C° (the space of compactly supported 
function, C°-valued, of class C°°) with domain the Sobolev space H'^{R''-) (g) C° [19]. 
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First of all, let us notice that the formalism and all the general results about constant fiber direct integrals 
are still valid in our setting of matrix- valued operators. We refer to \20\ Section XIII. 16] for a complete 
presentation of these results. 

For y € M*^, we denote by Ty the operator of translation by y which is defined, for u G L^(]R'^) ® and 
X € M'^, by {Tyu){x) = u{x — y). Then, because W is Z'^-periodic, the operator H is invariant by conjugation 
by r„, for every n G Z'^: 

Vn € Z"*, TnOHoT*=TnOHo T_„ = H. 

Thus, H is a Z'^-periodic operator. Let (ei, . . . , e^) be the canonical basis of W^. We recall that Co can be 
considered as the fundamental cell of the lattice Z'^, 

IfC^ is the fundamental cell of the dual lattice (Z^)* ~ 2ttZ'^, then is identified to the torus T* = W^/27rZ'^. 
Let e eT*. We denote by Vg the space of C^-valued, 0-quasiperiodic distributions in W^, which is the space 
of distributions u e P'(M^) (g) C° such that, for any n £ Z'^, T„?i = e-'"-^u. Let Tie = (LlcO^"^) ® C°) n Pg, 
endowed with the norm on L^{Co) (E) C°. We also define, for k e Z, the spaces = (H^^^iR"^) (E) C°) n V'g, 
where Hl'^ciW^) ® C° is the space of distributions that locally belong to the Sobolev space 
order to define the Fourier decomposition we will use later, it remains to define the space: 

■H = {ue (Lfo^(]R°') (g)L2(T*)) I y{x,e,n) G M'^ x T* x Z'^, u{x + n,e) = e'''-^u{x,e)} , 



Co = < xid + . . . + XdCd 



endowed with the norm: 



Vn G n, \\u\\n = ZZTT^ I M-^^)\\h{Co)^C^ 



vol(T*) 

For 6* G M'^ and n G L'^{R'^) O C°, we define Uu eUhy 

(2.13) Vx G R'^, ye G T*, {Uu)ix, 0) = ^ e^"-'^(r„n)(x) = ^ e'"-%(x - n). 

Actually, the expression (|2T3]1 is well-defined for u G S(R'^) (g) C°, the Schwartz space, and by Parseval 
theorem, this expression can be extended as an isometry from L^(M'^) (g to Ti. For every v ^Ti, we can 
define U*, the inverse of U by: 

(2.14) Vx G {U*v){x) = j^^ v{x, 9) d9. 
Indeed, we have, for v £ Ti and x G R'^, 

{UU*v){x) = J^e--^(C/*t;)(x-n)= J]e--^^^^j^ / v{x-n,9)d9 

= E "^^''ijn^) L ^'"■'^(^'^) = E e'"-'*n(x) = V{X,9). 
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Thus, U* is a left inverse for U which is an isometry from L^(]R"') C'^ to therefore U is unitary and 
U* = U"^. To obtain a Floquet decomposition for the operator H, it remains to prove that the operators 
U and H commute. As, for any j G {1, . . . , D} and any n G Z'^, the partial derivation dj commute with the 
translation r„, we have [dj, U] = 0. Thus, [— (g) Id, U] = 0. Then, using the Z'^-periodicity of W, we also 
have, for every u G L^(M'^) ® C°, every x G M'^ and every 6* G T*, 

{U oW){u){x,e) = ^ e^"- V(x - n)u(x - n) = ^ e^"-^VF(x)u(x - n) 

= VF(x) e^"-%(x - n) = VF(x)(?7n)(x, 0) = (VF o U){u){x, 6), 

as at X fixed, the multiplication by l^(x) G Sd(M) is continuous. Thus, [VF, C/] = and we finally have 
[i/, \J\ = 0. As we can see, even in the matrix-valued case we still have that H and U commute. Following 
[20j . we deduce that H admits the Floquet decomposition: 

r® 

(2.15) UHU* = Hg de, 

JT* 

where Hg is the selfadjoint operator H acting on T-Lg with domain Tig. Having this Floquet decomposition, we 
can continue to follow [2D] to obtain that Hg has a compact resolvant. It is a consequence on the assumptions 
made on the L^-regularity of W which ensure that H is elliptic. As Hg has a compact resolvent, its spectrum 
is discrete and we denote by 

Eo{e)<Ei{e) <...<Ej{e) < ... 

its eigenvalues, called the Floquet eigenvalues of H. Moreover, the functions 9 i— )• Ej{6), for j G N, are 
continuous and, if j tends to infinity, then Ej{6) tends to +oo, uniformly in 9. Actually, as Hg depends 
analytically on 9, we also have that 9 i-^ Ej{9) is an analytic function in the neighborhood of any point 
9^ G T* such that Ej{9^) is an eigenvalue of multiplicity one of ^^^o. 

The set Ej{T*) is a closed interval called the j-th spectral band of H and the spectrum of H is given by 

a{H) = [_}E,{T*). 

j&i 

If d > 2, the bands can overlap, but it is not the case in dimension 1 (except maybe at an edge point). 

Remark 2.1. Heuristically, the difference between the usual scalar-valued case (D = 1) and the matrix- 
valued case is that they are "D times" more Floquet eigenvalues in the matrix-valued case and thus the 
multiplicities of the Ej (9) 's are a priori bigger in the matrix-valued case than in the scalar-valued case. 



We finish this section by proving a result of non-degeneracy of the minimum of the first Floquet eigenvalue. 
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Proposition 2.2. Let 9^ G T* be a minimum of 9 i— )• Eq{6). Let v^{-,9^) be a Floquet eigenvector associated 
to the Floquet eigenvalue Eq[9^). Then, there exist 6 > and C > such that 

(2.16) ye € T*, \9-9^\<6 =^ C\9 - 9^\^ < Eo{9) - Eo{9^) < ^ ■ \9 - 0°^. 

This proposition means that 9^ is a nondegenerate minimum of the function 9 ^ Eq{9). 

Proof. We win follow the ideas of [8] and adapt them to the matrix- valued case. Let S^{-,9^) € GLd(C) be 
the fundamental solution the differential system 

(2.17) HsoS^{;9^) = Eo{9^)S^{;9^). 

Then, we define a new scalar product on the space Hgo by: 

yf,g G -Heo, < f,g >eo= [ < S''{x,9^)f{x),S''{x,9'')g{x) >cd dx. 

J Co 

We denote by T-lgo the Hilbert space (I-Lqo, < •, • >eo). Then, for every 6* E T*, we define on T-Lgo the operator 
Hg of domain 

D{He) = {u£neo \ S"^ {■ , 9°)u £ Uj} , 

and given by 

V^eT*, yueD{He), VxeR'^, {Heu){x) = S°{x,9°y\He - Eo{9^)){S°{x,9^)u){x). 

We have, using ()2.17p and integrating by parts, 

y9€T*,yu€D{Hg), <u,Hgu>eo = [ < {He - Eq{9^)){S°{x,9^)u){x),{S^{x,9^)u){x) y^i^ dx 

Jco 

< (5°(x,OVn)(x),(5°(x,0°)Vu)(x) >cd dx 

(2.18) = ||Vu||^o. 

Using general results on Dirichlet forms from [5l Appendic C], we have 

(2.19) Eoi9) - Eoie^) = inf 1^^^ ] <u,l >go= o] . 

I \MgO J 

Let u{x) = (0, . . . , 0, e'(''-^°)^, . . . ,0) e ngo. For this particular u, 

llVnll^o = ||(^ - Oe^(^-^")ll^o = |0 - e^PlHI^o, 
||n||go = 1 and < u,l >0O= 0. So, setting this u in (j2.18p and using (j2.19p . we obtain 

2 



Eo{9)-Eoi9')<ll^<\9-9r. 



For the lower bound in (|2.16p . we set: 



m+ = max ||5°(x, 0'^)|| and m- = min \\S^{x,9^) 
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Let u be such that ||u||go = 1 and < u,HqU >go being minimal where Hq is the operator constructed the 
same way as Hq from Hq, but in the case were the periodic potential W is equal to 0. In this case, Hq acts 
as — on D{Hg) = T-Lqo and, by Taylor formula, there exists 6 > such that, for \9 — < 6, 

< H^u, u > 
, \u\ 

Moreover, for every u € 'H'^ 



610, 

<Mu,u> < Vii, Vii >cD dx ^l|Vu||^o m\\\Vu 



H? /(^^^ < n, n >cD dx ^^^Iklleo H^Heo V"^- 

So we obtain the lower bound for C = ( ) . □ 



2.2. Wannier basis. We recall concepts used in [TOl E]. Let £ C L^(M'^) (g) be a closed subspace, 
invariant by Z'^-translations, i.e., for every n G Z'^, = r*n^Tn, where is the orthogonal projection on 
£. 

As is Z'^-periodic, it admits a Floquet decomposition similar to the one of H and, using the orthogonality, 
one gets: 

I-® 

u^= nf de, 

Jt* 

where Hg is the operator acting on Tig. The operator Hg is therefore an orthogonal projection acting on 
L'^{Cq) ® C^. As for {Hg)g^j*, the family {IIq)^^^* is continuous in 9 and thus is of constant rank. If we fix 
9 G T*, we can find an orthonormal system (wm,o)rrteAf) with M C N a set of indices independent of 9, that 
spans the range of 11^. Taking the image by U* of this orthonormal system, one gets an orthonormal system 
{■Wm,o)meM- If we set, for n G Z^, Wm,n = rn{wm,o), then {■Wm,n){m,n)eMxZ'i is an orthonormal basis of £. 
Such a system is called a Wannier basis of £. The vectors (5'm,o)mgM are called the Wannier generators of 
£. 

Let £ C L^(M'^) (g) be a space which is invariant by Z'^-translations. The closed subspace £ is said to be 
of finite energy for H if U^HU^ is a bounded operator. In this case, £ admits a finished set of Wannier 
generators. We now assume that £ is of finite energy for H. 

Let Jo be the set of indices of the Floquet eigenvalues of H which take the value for some values of G T*. 
We identify Jq to {1, ...,no}. Let Z be the set of 6* G T* for which there exists j G Jo such that, Ej{9) = 0. 
When 9^ is a nondegenerate minimum of Ej, Z is a set of isolated points (see |10j). It occurs when the density 
of states n has a nondegenerate behavior at (see [9]). For j G Jq, we define Zj = {9 £ T* ; Ej{9) = 0}. 
The sequence {Zj)j^jg is decreasing for the inclusion and Zi = Z. For 9^ G Z, Mqo C N is the set of indices 
such that Ej{9^) = 0. 



LIFSHITZ TAILS FOR MATRIX- VALUED ANDERSON MODELS 11 

We will denote by Wj{-,9) a Floquet eigenvector associated with the Floquet eigenvalue Ej{6) of H. For 
(^,^0 G (T*)2, we define Tg^e' ■ ^ He' by: 

Lemma 2.3. There exists {vj{-,9))j^jg, a family of functions on Tig, such that: 

1 ) for 9^ ^ Z and j G Mgo , there exists Vgo a neighborhood of Vqo in T* such that the map 9 G Vqo i— )• 
Vj{-,9) G Hg is real analytic (i.e. 9 i— Tg^goVj{-,9) is analytic as a function from Vgo to Hgo) and, for 
9 G Vgo, s^?m{{vj{-,9))j(,Mgo) = span((ii;j(-, 6l))jgAfgo)- 

2) For 9 G T*, the system {vj{-,9))ji=Mgo is orthonormal inTig andspan{{wj{-,9))j(zjQ) = span{{vj{-,9))j(zjQ). 
Proof We refer to [inilll], Lemma 3.1. □ 

In the next section, we will use this notion of Wannier basis and the notations we have just introduce to 
reduce our problem on estimating N(E) — A^(O^) to a discrete problem. 

3. Reduction of the problem 

The goal of this section is to give an estimate of N(E) — A^(O^) for an energy E close to 0. This will 
be accomplished by means of the IDS of certain reference operators, which are discrete operators. In this 
section we will use the notations introduced in Section [2j 

3.1. Reduction to a discrete problem. The reduction procedure consists into decomposing the operator 
according to various translation-invariant subspaces. The random operators thus obtained are what we 
consider as reference operators. They will be used to prove the upper bound on the IDS. 
We denote by ILo{9) the orthogonal projection in Tig on the vector space generated by {wj{-,9))j^jf). One 
defines 

(3.20) Uo = U-^(^j Uo{9)d9^U : ^^(M'^) ® C° ^ l2(M^) C°. 

Ho is an orthogonal projection on L^(M^) (g) and, for every n G Z,'^, we have r*noT„ = IIq. Thus, IIq is 
Z^-periodic. We set Sq = no(L2(M^) ® C°). This space is translation-invariant because of the Z'^-periodicity 
of IIq. Moreover £q is of finished energies for H as defined in p.6p . The main result justifying this reduction 
procedure is the following theorem which compares E i— )• N(E), the IDS of Hi^, to E ^£o^ the IDS of the 
discretize operator = UqH^Hq. 

Theorem 3.1. Let he defined by lil.5\) with the assumptions (HI), {H2) and {H3). There exist e > 
and C > 1 such that, for < E < e we have 

(3.21) < N{E) - N{0+) < Ng^iC ■ E), 
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where N^^ is the IDS of the discretize operator = HqH^Hq. 

Proof. See Theorem 4.1 in |10] . □ 

3.2. Periodic approximations. In order to get bounds on the density of states of no-f^cjllo, we wih now 
define periodic approximations of the operator H^i- For these approximations, we wih be able to control the 
density of state near by comparing it to some reduce operators. Then by taking a limit on the density 
of state of the reduced operators, we can get bounds on the density of states of HqHi^Uq and thus on the 
density of states of itself by using Theorem 13.11 Let /c G N*. We define the following periodic operator 

/ J"Wiix-in+l3)) 

(3.22) H^^k = -Ad®lD + W{x)+ Yl ^ \ 

neCfenZ-^ /3G(2fc+l)Z \ , 

The operator H^ f^ is (2/c + l)Z'^-periodic and essentially selfadjoint. It is an ff-bound perturbation of H with 
relative bound zero. Because of the {2k + l)Z'^-periodicity, we introduce the torus T^, = M'^/(2(2A; + l)7rZ'^). 
We also define A^^^ k, the IDS of k by 



de. 

where E^^kj is the j-th Floquet eigenvalue of the periodic operator H^^k- Let dN^^k be the derivative of 
A'^^fc, in the distribution sense. Ks E ^ Nuj,k{E) is an increasing function, dN^^k is a positive measure, it is 
the density of states of Hi^j^. Then, by [IHIISQ], for every ip £ C^(M), the distribution dN^^^k verifies 

(3.24) (<^,diV ) = ^/" tvnMH^Ae))d0,= ^^tr(lcMHu.,knc, 

[2tt}"- Jg^j* VOl(Cyt) V 

where tr(^) is the trace of a trace-class operator A. We index this trace by Tig if the trace is taken in Tig 
and here, the operator lc^.y:>{H^,k) is a trace-class operator. The proof of (|3.24p is given in [lOl Proposition 
5.1]. 

We want to take a limit on the density of states diV^ ^ of the periodic approximations in order to recover 
properties of the density of states of from properties of dN^^},. The following theorem ensure that it is 
possible. 

Theorem 3.2. 1) For any (p € Co°(M) and for almost every uj ^Vt, we have 

lim {ip,dN^^k) = {^,dN). 

k—^oo 

2) For any E & M a continuity point for N , we have lim E,{Ni^^kiE)) = N{E). 

Proof. The result of Theorem 13.21 is close to that of Theorem 5.1 of [lOj. The proof is also similar and is 
based on functional analysis. □ 
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4. Proof of theorem 11.51 

We will proceed in two steps. First, we will prove a lower bound and then an upper bound. 
4.1. Lower bound. In this subsection we prove 

Theorem 4.1. Let H^^, he the operator defined by /il.5\) with the assumptions [HI), {H2) and {H3). We 
have 



log 



(4.25) liminf 



log (n{E) - N{0+) 



d 

> 



E^o+ log E 2 

The proof of Theorem I4.lt As is the bottom of the spectrum, for e > we have N{e) — N{0) = 
N(e) — N(—e). To prove Theorem 14. H we will lower bound N(e) — N(—e). Then, for L large, we will show 
that H^^Cl recall that Hu,Cl restricted to with Dirichlet boundary conditions) has a large 

number of eigenvalues in [— e, e] with a large probability. To do this we will construct a family of approximate 
eigenvectors associated to approximate eigenvalues of H^^Cl ™ These functions will be constructed 

from an eigenvector of — A^C^/d + M^(x) associated to 0. Locating this eigenvector in 9 and imposing to 
to be small for n in some well chosen box, one obtains an approximate eigenfunction of H^^Ql ■ Locating the 
eigenfunction in x in several disjointed places, we get several eigenfunctions two by two orthogonal. 
In order to simplify the notations, we assume in what follows that = is a point where Eq{9) reaches 
0. From the same arguments as in [H] and using Proposition \2.2\ there exists C > such that, for 
/(.,6I) := (/i(-,6'), • • • ,/z)(-,6')) = vi{-,9) in L^{W^) (g)C°, {vi is the vector constructed in LemmaES]) one 
has 

(4.26) \\i-Ad^Ij, + W{x))fi;9)\\L2^Co)<^c^<C\9\^- 
We assume, without loss of generality, that /i ^ and we set 

(4.27) /(.,^) = Zl^(l,0,... ,0). 

|t7l| 

Let < ^ < 1 be a small constant. Let x £ C'o°(M) beeing positive, supported in [|,^] and such that 

X{tfdt = 2. 
For e > 0, we define 

d 

(4.28) We{9) = e-'"*\{x{e-"^e,)eL\T) and Wl{-,9) = We{9) ■ f{;9) L\Co) ^C"" . 

i=i 

Now let us estimate j|(-Arf (g) Id + VF(x))>V/[||^. We have 

||(-A, Id + TV(x))W/||2, = ^^^^ [ I|(-A, ® Id + W{x)){9)f{; 9)\\l,^^Co)^^o\We{9)\''d9. 



T* 



14 



H. BOUMAZA AND H. NAJAR 



So using (fOe]) and (g^H]), we get 

(4.29) ||(-Arf®/D + Ty(x))W/l|l^ < / \9\^\We{e)fde<C^e^ [ \e\^T\x\Oj)d0 

if ^ is small enoug h. For (3 G Z'^, we define 

^) = e-^'^-^W/l.e) and Wi,_^^^(-,0) = e-'^-^(n^^(^)W/)(-, 0), 
where Aci(C) is the cube defined by 

A„(C) = {n G z 
and nA^((^) is the orthogonal projection on Aa{C). 



^8' 



for l<j<d, \nj\ < 



We set 



Ulp{x)= j Wlp{x,e)de and = / 

JT* JT* 



For L large and /3 and (wi"^)„gzd well chosen, Uaepc, ^^^^ approximate eigenfunction of H^^^Cl 

associated to an approximate eigenvalue in the interval [—£, e] ■ 



if 



We notice that U. 



f 



G L 



2/Tn>d\ 



' C° and Ul]^^^^^ G L'^{R'^). As in [14j one gets that 



f 



> m 



fi 



Now we have to look to the conditions under which we have 



> C7 > 0. 



Arf ® Id + W{x) ]uj 



L2(Rd),g,CD 



Note that 



(4.30) 



f 



L2 



< 



< 2 



+ 2 



Equations (j4.29p give the bound on the first member of (|4.3ip . It just remains to control the second term. 
To do so, one needs the following lemma 



Lemma 4.2. Let ( = £2 . There exists K > 0, such that 

2 



(4.31) 



f 



L2 



<e'' + K\ A„(e) 



(n) 

sup OJl 

ne/3+2Atj{e) 



Before proving this lemma let us use it to finish the proof of Theorem 14.11 
Taking ()4.29p and (|4.3ip into account, we get that there exists K > such that 



(4.32) 



f 



2 



< — + K[ sup ml 

ne/3+2Ac,(e) 



(n) 
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Now, for L large, we may divide Cl into disjoints cubes of size 2Ka{e). For a < ^, there exists (7 > 
such that L{e) satisfies 

We can find /3i, . . . , /^^(e) in lY such that : 

|J(/3,+A„(e)) cCi and, for jV/, (/3,- + 2A„(£)) n (/3,v + 2A„(e)) = 0. 
i=i 

It particular, for j ^ j' , U^^ ^ and ^/^^ ^ are orthogonal. Then, 

e(# I eigenvalues of Hc^H^Hcl in > e {1, • • • | ||i^.;Z^i,,^^,,J|L2(iRd)^CD < ^}) 

(Lie) \ 

(4.34) > e( J^i?,(a;)j , 

where 



1 if K I Aa(e) |2 • (sup„g^^.+2A<,(e) '^i"^) < T- 



if not. 

The iBj)i^j^x,{e) i.i.c?. Bernoulli random variables. So equations (j4.34p and (j4.33p imply that there 
exists C > such that one has 

^2L + iy K*{ ^ig^^^al^^s oiUc,H^Uc, in [-e,6]}) > ^^^^P{B, = 1) > ^6'^P(Si = 1). 
Hence, taking the limit L — > oo, we get that, for e > small, 



(4.35) N{e) - N{-e) > ^e''P(i?i = 1). 



(n) . e • e 



It just remains to estimate P{Bi = 1). If, for 1 < j < L{e) and n € /3j+Aa{s), one has uj\ < / ^ 12 ' 

zA I i\q,[E) I 

then for e rather small 

i^|A,(e)|2( sup JY^y<^^ 
As the random variables are i.i.d., one has the estimate 



2 



Hence, taking the double logarithm of (j4.35p . using assumption (H3) and the fact that j^ka(e) = £"(2+'^")^ 
we get that 

log log (iV(e) - Af(0)) ^ 

(4.36) lim ^ i— > da. 

^ ' 6^0+ loge - 2 

The equation (I4.36P is true for any a > 0, by letting a tend to 0, we end the proof of Theorem 14. 1[ □ 
It remains to prove Lemma 14.21 to finish this section on the lower bound. 
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The proof of Lemma 14.21 We have : 



(4.37) 

Then, 
(4.38) 



f 



L2(Co)(g>CD 



Here we used the fact that lA^^ „ , /, and U^^a are close to each other. 



dx. 



We set 



S, 



13,6 



dx. 



As Vi is supported in Co, one gets that there exists K > {) such that 

2 



(4.39) 



dx. 



One sphts the sum on rj (^U^ in two parts according to whether r] belongs to /3 + AQ,(e) or not. 
• For the sum on r/ ^ /3 + (e) , there exists C > such that 



(4-40) E E -1 

r?^^+A„(_) \nGr?+A„(£) 



(n) 



Co 



e,ryv 



dx. 



By the use of the non-stationary phase and following the same computation done in [T3], one proves that, 
for m an integer, there exists > such that 



(4.41) 



dx] < K^e'^"'. 



• For the sum on r/ G /? + Aq, (e) , one gets 
(4.42) 



E E 

And there exists -fC > such that 



Co 



z^i;i,(x) 



dx < |AQ,(e)| • sup u}\ 

\ 7ie/3+2Ac(e) 



dx. 



(4.43) 



dx < K. 



We finally use (09]) . (lOOD . (OT]l . ([02]l and (03]) to get that 



<e^ + K\ AJe) 



(n) 

sup UJ\ 
nG/3+2Ac(e) 



So we obtain Lemma l4.2i D 

To finish the proof of Theorem 11.5] it remains to prove the upper bound. 
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4.2. Upper bound. We start this section by recahing that, as we deal with the bottom of the spectrum, we 
have non-degeneracy of the first Floquet eigenvalue at the bottom of the spectrum as shown in Proposition 
12.21 Using this, we prove the following theorem: 

Theorem 4.3. Let he the operator defined by il.5\) with the assumptions (HI), (H2) and (H3). Then 

\og\\og{N{E)-N{Q+))\ d 
limsup ; < . 

E^o+ ^ogE - 2 

The proof of Theorem 14. 3L To prove the upper bound, it is enough to prove the same upper bound on 
Ngf^ (as defined in Theorem 13. ip . To do this, we show that Ng^^ (and so N) may be compared to the IDS of 
some well chosen discrete Anderson model whose behavior of its IDS is already known. 
We begin by isolating the contributions from the various points for which Ej{6) take the value 0. We recall 
that the band at is generated by {Ej{9))i<j<no- For 1 < j < no, Zj = {6 ^ T*;Ej{9) = 0}. The sequence 
(Zj)i<j<„Q is decreasing (Zj^i C Zj). Let 9^ £ Z. We set j{9^) = supMgo with Mgo = {j ; 1 < i < 
no,Ej{9^) = 0}. We replace the Floquet eigenvectors (•, 6'))i<j<j(go) associated to (£'j(^))i<j<j(go) by 
the vectors ^))i<j<j(6»") constructed in Lemma [2T3l They are analytic in a neighborhood Vqo of 9^. Let 
9 be close to 9^. The operator H^{9) = no-ff(^)no is unitarily equivalent to the multiplication operator by a 
function on L^(T*) with values in A^„g(C). This matrix-valued function takes the following block diagonal 
form : 

fBj(^gO){9) ... \ 

Ej^eo)+i{9) ... 

■■. ■•. : 

V ... s„oW/ 

where the matrix Bj(^Q0-^{9) is of size j{9^) x j{9^) and is given by 

/ {vi{-,9),H{9)vi{-,9))l2(^Co)(&c° ■■■ {M-,^), Hi9)vj^go){-,9)) l2(^Co)^c° 

\{vj(eo)i-,0),H{9)vi{-,9))i^2(^Co)(E)C^ ... {vj(^eo){-,9),H{9)vj^go)i-,9))L2(^Co)(^c^, 
The matrix Bj(^0o-^{9) has (£'j(6'))x<j<j(6iO) for eigenvalues. The operator = IIoV^IIo is unitarily equivalent 
to the multiplication operator by the matrix with entries ((V^^j, Wj))i<jj<no- 
For u G L2(T*) (g) CD, Hon = Y^Zii'^'^'^i) L^m^^c^-^i- For 9^ G Z, we set 

d 

^eo(^) = 5](l-cos(0,-e°)). 

i=i 

We recall that the eigenvalues (£'j(0))i<j<j(5io) are non-degenerate at 0. So there exists Vgo (an open 
neighborhood of 9^) and C > 1 such that, for 9 G Vgo, we have, for 1 < j < ji9'^), CEj{9) > Wqo{9) and, 
for j > j{9^), CEj{9) > 2. We remark that the neighborhood Vgo can be chosen such that Vqo C Vqo, where 
Vgo was defined in Lemma [2^ 
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Let Hgo{9) be the uq x no diagonal matrix with identical diagonal entries equal to zugo. For 9 G Vgo, we have 

(4.44) H^goiO) <C ■ H^e). 

Finally, we note that {Vgo)0O(zz can be chosen so that they cover T*, (i.e. Ugog^Vgo = T*) and such that 
each one of them contains only one point of Z (i.e. for 6 ^ Z,6' ^ Z such that 9 ^ 9', we have 0' ^ Vg). We 
order the points in Z = {9^; 1 < k < mo}, where niQ = i^Z. Let (xfc)i<fc<mo be functions in C^(T*) which 
form a partition of the unity on T* such that, for I < k < uiq, supp (xk) C Vgk, < Xfc < 1 and Xfc = 1 in a 
neighborhood of 9^. 

So there exists C > 1 such that, for any € T*, we have, 

, mo mo 

(4.45) —<Y^^l<i and ^ H'^,{9)xI < CH^i9). 

k=l k=l 

For t G (L^(T*)(g)C°)(XiC""(X>C™", we set t = {tj^k)i<j<no:i<k<mo- We consider t as a system of mo columns 
denoted by (t.,fc)i<fc<mo- Each column belongs to (L2(T*)(8)C^)(g)C"". We endow (L2(T*)(g)C°)®C"«(g)C™o 
with the scalar product generating the following Euclidean norm: 



'"On 2 

EWt-k = 



(L2 (T* ) iX)CD ) iX)C"o ® C™o 



(L2(T*)®CD)®C"o 



tj,k 



k=l ^^j^'i^O:^^k<mo 

We define the mapping S : (^^(T*) C°) (g) C"" — > {L^{T*) C°) (g) C"« C^o by 



2 



L2(T*)®CD 



5(t) = {Xkt)l<k<mo = iXktj)l<j<no, l<fc<mo, if t = (*i)l<j<n.o ^ (L2(T*) (g> C°) «) C"«. 

Here, for any 1 < j < no, tj = {tij)i<i<D G L'^{T*) O C°. 

The adjoint of 5, S* : {L'^{T*) C°) C"« C"^" ^ (^^(T*) O C°) ® C"o is defined by 

S*{t) = I Yl ^khk) for t = (t,-,fc)i<,-<„o;i<fe<mo e (^'(T*) ^ C°) ® C"« ® C"«. 

V l<A,<mo / l<j<no 

Here, for any 1 < j < no and any 1 < k < niQ, we have tj^^ = {tij^k)i<i<D- According to equation (|4.45p 
we have ^/ < 5* o 5" < /, (here / is the identity in {L^{T*) (g) C°) (g) C""), thus 5 is one to one. Using the 

(n) 

boundedness assumption on the Vi and on the support of the w^- , one shows the following lemma: 

Lemma 4.4. There exists C > such that, for t G {L^{T*) (g) C°) (g C"^, we have 

(-^(2'S'(0) 'S'(i))(L2(']r.)0CD)(g)C"OiX)C'"o < (^(-f^i!^) *)(L2(T*)iX)CD)iX)C"o ) 
where the operator acting on (L^(T*) (g) C^) (g C"" (g) C"*" zs defined by 

-^(2* = ( Hlti j k + Vi^^Ati j k I 

V ' /l<i<-D, l<i<no, l<fc<mo 

iJere, is the multiplication by Wgk acting as a multiplication operator on L^(T*), V^^ = YlneZ'' ^i^^^n, 
where n„ is the orthogonal projection on the vector 9 ^ e'"^ in L^(T*), and is defined in Section [3A[ 
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The proof of this lemma fohow the same steps as Lemma 5.5 in [14J. We use it to end the proof of Theorem 
14.31 Lets first notice that the operator could be written as a direct sum of no copies of mg x D random 
scalar-valued continuous Anderson models. Indeed, we can write 

(l2(t*) (g)C°) (»c"« (g)C™" = L^{T*)(g)a0O 

l<i<D, l<i<no, l<fc<mo 

Here, for 1 < j < I, we use the notation = {0}-'^^ x C x {0}'"-' . So H^^ is unitarily equivalent to 

(8) 

l<i<D, l<j<no, l<fc<mo 

Here H^^'^ acts on L^(T*) (8) C* (E) O (S> C^. Using the discrete Fourier transformation, we get that for every 
A; € {1, . . . jmo}, H^^'j, is unitarily equivalent to h^^'^, where h^'^'^ acts on I'^iWj'^) and is defined by 

(4.46) <f = -A^. + ^ 

Here, if 5n is the vector {5'^)p^id where (5,^ is the Kronecker's symbol, then 7r„ is the orthogonal projection 
on 5n and — A^d is the discrete Laplacian defined by : 

(4.47) Vn G l\'L''), (A^,«)„ = ^ ^ {un- u^). 

\m—n\=l 

Using the fact that for operators A and B, we have A^(^ ® B,E) = N{A,E) + N{B,E) (see ^l^)! we get 
that 

D 

(4.48) Ne,{e) < no x mo x ^ A^(<f , C.mQ.e). 



4 = 1 



To satisfy assumptions of [21], we set, for every i G {1, . . . , D}, Sj = sup„g^d w^-"''' and : 

1 S./2 it 4 £(»./2.0.1 

By changing ujf^^ into tDj-"'' in (j4.46p . we define a new operator which we denote by h^'f^- We notice that 
^uH^ lower bound h^'^'^ with the same bottom of the spectrum. As it is known that each h^'^'^ exhibits 
Lifshitz tails with Lifshitz exponent —d/2 (see l6l[21j), using Theorem 13. II and (|4.48p . we get that 

log|log(jV(£)-jV(0+))| d 

limsup < . 

e^o+ log e 2 

This ends the proof of Theorem 14.31 □ 
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